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At HERA heavy quarks may contribute up to 30% of the structure function F2. The 
potential of including heavy-quark data in the ZEUS PDF fits is explored, using D* 
double differential cross-sections as well as the inclusive quantities ^2^^, -F^*^. The 
introduction of heavy quarks requires an extension of the DGLAP formalism. The 
effect of using different heavy flavour number schemes, and different approaches to the 
running of as, are compared. 

Parton Density Function (PDF) determinations are usually global fits [21 [31 [U [S] , which 
use inclusive cross-section data and structure function measurements from deep inelastic 
lepton hadron scattering (DIS) data as well as some other exclusive cross-ections. The kine- 
matics of lepton hadron scattering is described in terms of the variables Q^, the invariant 
mass of the exchanged vector boson, Bjorken x, the fraction of the momentum of the incom- 
ing nucleon taken by the struck quark (in the quark-parton model), and y which measures 
the energy transfer between the lepton and hadron systems. The differential cross-section 
for the neutral current (NC) process is given in terms of the structure functions by 

= ^ [Y+ F,{x, Q') - F,{x, Q^) T xF,{x, Q')] , 

where y± = 1 ± (1 — j/)^. In the HERA kinematic range there is a sizeable contribution to 
the F2 structure function from heavy quarks, particularly charm. Thus heavy quarks must 
be properly treated in the fomalism. Furthermore fitting data on charm production may 
help to give constraints on the gluon PDF at low-x. 

The most frequent approaches to the inclusion of heavy quarks within the conventional 
framework of QCD evolution using the DGLAP equations areS 

• ZM-VFN (zero-mass variable flavour number schemes) in which the charm parton 
density c{x,Q^) satisfies c{x,Q'^) — for < and nf = 3 + 0{Q^ — fj,^) in the 
splitting functions and /3 function. The threshold /i^, which is in the range rri^ < 
/Lt^ < 4m^, is chosen so that F|(a;,(5^) = 2e'^xc{x,Q^) gives a satisfactory description 
of the data. The advantage of this approach is that the simplicity of the massless 
DGLAP equations is retained. The disadvantage is that the physical threshold = 
Q^ij — 1) > 4m^ is not treated correctly {W is the 7*3 CM energy). 

• FFN (fixed flavour number schemes) in which there is no charm parton density and 
all charmed quarks are generated by the BGF process. The advantage of the FFNS 
scheme is that the threshold region is correctly handled, but the disadvantge is that 
large ln((5^/m^) terms appear and charm has to be treated ab initio in each hard 
process. 

• GM-VFN (general mass variable flavour number schemes), which aim to treat the 
threshold correctly and absorb ln((5^/m^) terms into a charm parton density at large 
Q^. There are differing versions of such schemes [6l[7] 



^Charm production is described here but a similar formalism describes beauty production 
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Figure 1: Comparison of predictions for Jlf'^, from fits which use the GM-VFN scheme 
and the PEN scheme with two different factorisation scales: on the left hand side the FFN 
schemes still use a VFN treatment of a^, whereas on the right hand side a 3-flavour as is 
used. 



For the ZEUS PDF analyses [HI [5] , the heavy quark production scheme used was the 
general mass variable flavour number scheme of Roberts and Thorne (TR-VFN) [HI [7]. 
However we also investigated the use of the FFN for 3-flavours with the renormalisation 
and factorisation scale for light quarks both set to but the factorisation scale for heavy 
quarks set to + 4™^. The reason for these choices of scheme and scale is that these are 
the choices made in the programme HVQDIS [10] which was used to extract Fl'^ from data 
on D* production. Furthermore, it has recently become evident that the use of the FFN 
scheme implies a treatment of the running of ots which is different from that of the VFN 
schemes. In VFN schemes as is matched at flavour thresholds [T^], but the slope of as is 
discontinuous. In the FFN scheme we must use a 3-flavour as which is continuous in . 
This requires an equivalent value of as[Mz) = 0.105 in order to be consistent, at low Q^, 
with the results of using a value of as[Mz) = 0.118 in the usual VFN schemes. 

In Fig [1] we compare different heavy quark factorisation scales and different treatments 
of the running of as for predictions of i^^'^- We see that within the FFN scheme the choice of 
the heavy quark factorisation scale makes only a small difference at low Q^. The treatment 
of as gives larger differences. The FFN scheme and GM-VFN scheme differ for almost all 

if as runs as for the VFN schemes. However if a 3-flavour as is applied in the FFN 
schemes there is much better agreement of all schemes at higher . 

There is now new data on .F^f'^ [H] Fl** [15] from HERA-II running to add to older 
the HERA-I charm data [5]. To investigate the potential of these data to constrain the 
gluon PDF, we used the ZEUS-pol PDF fit [T3j and added the charm data. The new data 
do not influence the central values of the ZEUS-pol fit significantly. However there is a 
small improvement in the precision of the low- a; gluon. Fig. [2] compares the PDFs and their 
uncertainties, as extracted from the ZEUS-pol PDF fit, with the those extracted from a 
similar fit including the F^"^ and F2'' data. This illustrates that the charm data has the 
potential to constrain the gluon PDF uncertainties. 

We have also compared fits using the GM-VFN formalism with those using the FFN 
formalism with 3-fiavour as. When using the FFN formalism one should not really use 
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Figure 2: The u-valence, d-valence, Sea and gluon PDFs and their fractional uncertainties 
at Q2 = lOGeV^, from a) the ZEUS-pol PDF fit (left) and b) a smilar fit with Ff and F^^ 
data included (right). □•production 
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Figure 3: Double differential cross-sections for D' prodmction. The red lines show the 
predictions of the ZEUS-S-13 NLO PDF fit using the Petersen fragmentation function for 
the Z?*, whereas the blue lines show these predictions using the Lund fragmentation function. 



high-Q^ data, because large ln(Q^/m^) terms are not resummed. Another short-coming is 
that the NLO FFN coefficient functions are not available for the CC processes. Since the 
CC reactions at HERA are at high-Q^ we use ZM-VFN coefficient functions. In practice the 

for these fits is not bad. The main difference between FFN and VFN fits is in sensitivity 
to the charm quark mass nic, with the FFN fits preferring a low value TOc = 1.35 GeV and 
the GM-VFN fits favouring a higher value rric ~ 1.45 GeV. 

The small impact of the heavy flavour data on the global fit may be because we are 
not using the charm data optimally. i^2^'^ is a quantity extracted from D* cross-sections by 
quite a large extrapolation. It would be better to fit to those cross-sections directly. The 
evaluation of the theoretical predictions involves running the NLO programme HVQDIS 
for each iteration of the fit. However, one can shorten this process by using the same 
method as was used for the ZEUS- JETS fit [11]. The PDF independent subprocess cross- 
sections are output onto a grid, such that they can simply be multiplied by the PDFs at 
each iteration. The data used are the nine double differential cross-section measurements of 
<Pa{D*)/dQ^dy [5], see Fig. [31 When fitting D* cross-sections, as opposed to and inclusive 
quantity like i^a^"^, one must use the FFN scheme since the prediction grids are calculated 
using HVQDIS. This means that we cannot use ZEUS high-Q^ data. Hence we chose to 
use the ZEUS-S global fit [5] as the basis for our fit, with a cut-off < SOOOGeV^. The 
parametrisation was slightly modified to free the mid- a; gluon parameter p^ig) and the low- a; 
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Figure 4: The gluon PDF and its fractional uncertainties for various bins Left: before 
D* cross-section data are input to the ZEUS-S-13 fit. Right: after D* cross-section data are 
input to the ZEUS-S-13 fit 

valence parameter P2(u) = P2{d), such that the paramctrization is like that of the ZEUS- 
JETS and ZEUS-pol fits. This fit is called ZEUS-S-13. Figure H shows the difference in the 
gluon PDF uncertainties, before and afer the D* cross-sections were input to the ZEUS-S-13 
global fit. Disappointingly the uncertainty on the gluon is NOT much improved. 

Should we have expected much improvement? There are two aspects of the fit which 
could be improved. The predictiond for the D* cross-sections have more uncertainties than 
just the PDF paramctrization. A further uncertainty is introduced in the choice of the 
c D* fragmentation The Petersen fragmentation function was used for the fit predictions. 
However, looking back at Fig [3] we can see that the Lund fragmentation function seems to 
describe the data better. To best exploit the charm data in future we need to address such 
aspects of our model uncertainty. Secondly, this study on the D* cross-sections used only 
the HERA-I charm data. We look forward to the 5-fold increase in statistics expected from 
HERA-II charm and beauty data. 
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